HISKP-TH-06/31 



Universal Properties of the Four-Body System with Large 

Scattering Length 

H.-W. HammeiB 

Helmholtz-Institut fur Strahlen- und Kernphysik (Theorie), 
Universitdt Bonn, 53115 Bonn, Germany 

L. PlatteiQ 

Department of Physics and Astronomy, 
Ohio University, Athens, OH 45701, USA 

(Dated: February 9, 2008) 

Abstract 

Few-body systems with large scattering length have universal properties that do not depend on 
the details of their interactions at short distances. We study the universal bound state properties 
of the four-boson system with large scattering length in an effective quantum mechanics approach. 
We compute the four-body binding energies using the Yakubovsky equations for positive and 
negative scattering length. Moreover, we study the correlation between three- and four-body 
energies and present a generalized Efimov plot for the four-body system. These results are useful 
for understanding the cluster structure of nuclei and for the creation of weakly-bound tetramers 
with cold atoms close to a Feshbach resonance. 

PACS numbers: 03.65. Ge, 36.40.-c, 21.45. +v 



* Electronic address: hammer@itkp.uni-bonn.de 



^Electronic address: lplatter@phy.ohiou.edu 



I. INTRODUCTION 



Effective theories are a powerful framework to exploit a separation of scales in a physical 
system. The long-distance degrees of freedom must be included dynamically in the effective 
theory, while short-distance physics enters only through the values of a few coupling con- 
stants, often called low-energy constants. Effective theories are widely used in many areas 
of physics. Recently, a considerable effort was devoted to applying effective field theories 
in nuclear and atomic physics [l|, H, Hj|. If there is no exchange of massless particles, any 
interaction will appear short ranged at sufficiently low energy. One can then use a very gen- 
eral effective theory with short-range interactions only to describe the universal low-energy 
properties of the system. Such a theory can be applied to a wide range of systems from 
nuclear and particle physics to atomic and molecular physics. 

In this paper, we focus on few-boson systems with a large two-body scattering length. 
They are characterized by an unnaturally large two-body scattering length a which is much 
larger than the typical low-energy length scale I given by the range of the interaction. Such 
systems display a number of interesting effects and universal properties that are independent 
of the details of the interaction at short distances of order I. The simplest one is the 
existence of a shallow two-body bound state (called dimer) with universal binding energy 
B2 = h 2 / (ma 2 ) + Oil I a) for positive a and where m is the mass of the particles. Low-energy 
observables can generally be described in a controlled expansion in l/\a\. In the two-boson 
system, the effective theory reproduces the effective range expansion but the structure of the 
three-boson system is much richer 0]. It has universal properties that include a geometric 
spectrum of three-body bound states (so-called Efimov trimers), log-periodic dependence 
of three-body observables on the scattering length, and a discrete scaling symmetry 0, lEf . 
These features can be understood as the consequence of a renormalization group limit cycle 
in the three-body system jl, 0, 0]- 

In the effective field theory (EFT) formulation of Bedaque et al. 0, 0], the limit cycle is 
manifest in the renormalization group behavior of a contact three-body interaction required 
for consistent renormalization. This implies that at leading order in l/\a\, the properties of 
the three-boson system with large scattering length are not determined by two-body data 
alone and one piece of three-body information (such as a trimer energy) is required as well. 
This three-body information can be specified in various ways jl, @, 0, H, @, [13, EH 12] • In 



the following, we will use the parameter L 3 introduced in [12| . As an alternative to the EFT 



formulation, one can construct an effective theory in a quantum mechanical framework [13 



The contact operators in the field theory are then replaced by an "effective potential" built 
from smeared out 5-function potentials and derivatives thereof. The qualitative features 
of the renormalization are very similar to the EFT formulation and also show the limit 
cycle behavior. In the case of positive scattering length, this approach has successfully been 
applied to the three- and four-boson systems p, 12, 14] and to the four-nucleon system (IBj . 
Here we extend this work to negative scattering lengths. 

The four-body problem has previously been studied in more traditional approaches. Early 
studies include the Schrodinger equation with separable potentials [13] and field-theoretical 
models with separable expansions of the three-body T-matrix [l8[. For a review of these 



and other early studies, see Refs. [19|, l20j. The nuclear four-body problem has recently been 
benchmarked by comparing various modern calculational approaches [2lJ. An overview of 
recent calculations for the four-body system of 4 He atoms can be found in Ref. (22| . 

None of these previous works focused on the universal properties of the four-body system 
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with large scattering length which remain poorly understood. A first step towards under- 



standing the four-body system with large positive a was taken in Ref. [12|. By means of an 
explicit calculation, it was demonstrated that the renormalization of the three-body system 
automatically ensures the renormalization of the four-body system in this case. Therefore no 
four-body parameter should enter at leading order. These results were applied to calculate 
the 4 He tetramer ground and excited state energies and good agreement with the Monte 
Carlo results of Blume and Greene was found 22|. Within the renormalized zero-range 
model, however, Yamashita et al. 23J recently found a strong sensitivity of the deepest four- 
body energy to a four-body subtraction constant in their equations. They motivated this 
observation from a general model-space reduction of a realistic two-body interaction close 
to a Feshbach resonance. The results of Ref. [l2| for the 4 He tetramer were also reproduced. 
Yamashita et al. concluded that a four-body parameter should generally enter at leading or- 
der. They argued that four-body systems of 4 He atoms and nucleons (where this sensitivity 



is absent [12|, |15|, [16J) are special because repulsive interactions strongly reduce the proba- 



bility to have four particles close together. However, the renormalization of the four-body 
problem was not explicitly verified in their calculation. Another drawback of their analysis 
is the focus on the deepest four-body state only. Therefore, it remains to be seen whether 
their findings are universal or are an artefact of their particular regularization scheme. 

The purpose of this paper is to extend the study of Ref. [l2| in two important ways. 
First, we carry out calculations for positive and negative scattering length and study the 
universal correlations between the three- and four-body binding energies. Second, we map 
out the scattering length dependence of the shallowest two four-body states and summarize 
the spectrum in a generalized Efimov-plot. We will follow Ref. 12] and work at leading order 
in l/\a\ using the framework of non-relativistic quantum mechanics to construct an effective 
interaction potential. The universal properties of the four-body spectrum will be useful 
for understanding the cluster structure of nuclei [24| as part of the planned experimental 
program at FAIR and for the creation of weakly-bound tetramers with cold atoms close to 
a Feshbach resonance. A first step in this direction was already taken in Ref. (2o| . 

The organization of the paper is as follows: In Sec. [Til we briefly describe our formalism 
for the two-, three-, and four-body bound state equations. The discussion of the universal 
correlations and the four-body spectrum follows in Sec. IIHI Finally, we summarize and 
present an outlook in Sec. IIV1 



II. FEW-BODY BOUND STATE EQUATIONS IN EFFECTIVE THEORY 

The effective low-energy interaction potential generated by a non-relativistic EFT with 
short-range interactions can be written down in a momentum expansion. In the two-body 
S-wave sector, it takes the general form 

(k'|V|k) = \ 2 + \ 2 ^k 2 + k' 2 )/2 + ... , (1) 

where k and k' are the relative three-momenta of the incoming and outgoing particles, 
respectively. Similar expressions can be derived for three- and higher-body interactions. 
The exact form of the potential depends on the specific regularization scheme used. The 
low-energy observables, however, are independent of the regularization scheme (up to higher 
order corrections) and one can choose a convenient scheme for practical calculations. 

We regularize the potential in Eq. (TjQ) by multiplying with a Gaussian regulator function, 
exp[— {k 2 + k' 2 )/A 2 ], with the cutoff parameter A. This factor strongly suppresses high- 
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momentum modes in the region k,k' > A where the effective potential is not valid. The 
cutoff dependence of the coefficients A2(A), A2,2(A), ... is determined by the requirement 
that low-energy observables are independent of A. 

The expansion in Eq. is only useful in conjunction with a power counting scheme that 
determines the relative importance of the various terms at low energy. The leading order is 
given by the A2 term which must be iterated to all orders, while the other terms give rise to 



higher-order corrections that can be included perturbatively [26|, |27 



In this paper, we will 

restrict ourselves to leading order in the expansion in l/\a\ and include only the A2 term. In 
the three-body system, a momentum-independent three-body interaction term A3 must be 
included together with A2 at leading order in l/\a\ [a M| . Effective range effects and other 
higher-order corrections can be included as well @,Tlo[ H 0, 0, H(J. In the four-body 
system, no new parameter enters at leading order in l/\a\ and only A 2 and A 3 contribute 

In order to set up our conventions and formalism, we will briefly review the bound 
state equations for the two-, three-, and four-body systems. For a more detailed discussion 
including explicit equations in momentum space, we refer the reader to Refs. [l2|, 31 . 



A. The Two-Body Sector 

We write the leading order two-body effective potential in momentum space as: 

(p|V|q> = (pb>A 2 (s|q> , (2) 

where A2 denotes the two-body coupling constant and q (p) are the relative three-momenta 
in the incoming (outgoing) channel. The regulator functions 

(p\g) = g{p) = exp(-p 2 /A 2 ) , (3) 

suppress the contribution from high momentum states. In the few-body literature, they are 
often called "form factors" . For convenience, we will work in units where Planck's constant 
H is set to unity: h = 1. 

The interaction ([2]) is separable and the Lippmann-Schwinger equation for the two-body 
problem can be solved analytically. The two-body t-matrix can be written as (32| : 

t(E) = \g)r(E)(g\ , (4) 

where E denotes the total energy. The two-body propagator r(E) is then given by 

(5) 

A two-body bound-state appears as a simple pole in the two-body propagator r at energy 
E = —B-i- Thus for a > 0, the two-body coupling constant \2{B2, A) can be fixed from the 
binding energy B 2 , which is directly related to the scattering length by a = \j \Ju1B2 up to 
higher order corrections in l/\a\. For negative scattering length, there is no shallow dimer 
state and the coupling constant A2 is fixed from matching to the effective range expansion 
for positive energies. 



t(E) 



1/A 2 -4tt / dqq 



mE — q 1 
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B. The Three-Body Sector 



The low-energy properties of the three-body system for a given effective potential can 
be obtained by solving the Faddeev equations. The full wave function can be decomposed 
into four components: one for each two-body subcluster and one for the three-body cluster 
[33|]. For identical bosons, the three-body wave function is fully symmetric under exchange 
of particles and the Faddeev equations simplify considerably. In this case, one only needs 
to solve equations involving one of the two-body Faddeev components and the three-body 
component. The two remaining two-body components can be obtained by permutations of 
particles. 



We follow Glockle and Meier [33j and decompose the full three-body wave function as 



$ = (l + P)^ + ^3, where P = P 13 P 23 + P l2 P 23 (6) 

is a permutation operator that generates the two not explicitly included Faddeev components 
from ip. The operator Pij simply permutes particles i and j. Since we are interested only in 
the binding energies and not in the full wave function, we can eliminate the component ip 3 
and obtain an equation for if) alone: 

ip = G tPi) + G tG t 3 (l + P)iP . (7) 

where Go denotes the free three-particle propagator. Moreover, t is the two-body t-matrix 
for the two-body subsystem described by the component if) and t 3 is a auxiliary t-matrix 
defined by the solution of the three-body Lippmann-Schwinger equation with the three-body 
effective interaction 

V 3 = |£>A 3 <£| , (8) 
only. The three-body regulator function |£) is defined as 

( Ul , u 2 |£) = exp (-(«? + 3^/4)/A 2 ) , (9) 

where Ui and U2 are the usual Jacobi momenta. Note that t 3 is only a technical construct 
that is generally cutoff dependent and not observable. For the derivation of an explicit 



representation of Eq. (j7]) in momentum space, we refer the reader to Refs. [12|, [31 . 

The three-body binding energies are given by those values of E for which Eq. ([7j) has a 
nontrivial solution. By expressing the two-body coupling constant A2 in terms of the binding 
energy of the shallow two-body bound state, we have already renormalized the two-body 
problem. The three-body problem can therefore be renormalized by requiring one of the 
three-body binding energies to be fixed as one varies the cutoff. All other low-energy three- 
body observables can then be predicted. This renormalization procedure determines the 
three-body coupling constant X 3 (B 3 ,A) uniquely. 

The dimensionless coupling constant A 4 As(A) shows a limit cycle behavior 12T] . For large 
values of the cutoff A, A A X 3 flows towards an ultraviolet limit cycle. For A — > 00, it has the 
limiting behavior 

wax c sin(g ln(A/L 3 ) - arctan(l/g )) 

A 4 sin(s ln(A/L 3 ) + arctan(l/s )) 1 ' 

where Sq ~ 1.00624 is a transcendental number that determines the period of the limit cycle. 
The constant c = 0.074 ± 0.003 was determined in Ref. The discrete scaling symmetry 



associated with a limit cycle is manifest in Eq. (fTUl) . If the cutoff A is multiplied by a factor 
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exp(n7r/so) ~ (22.7) n with n an integer, the three-body coupling A3 is unchanged. L 3 is 
a three-body parameter generated by dimensional transmutation. It can be determined by 
fixing a three-body binding energy P 3 . Alternatively, one could use a three-body binding 
energy directly to characterize the value of the three-body coupling A 3 at a given cutoff 

MM- 



C. The Four-Body Sector 

We now turn to the four-body sector. The four-body binding energies are given by the 
Yakubovsky equations which are a generalization of the Faddeev equations to the four-body 
system. The full four-body wave function \1/ is first decomposed into Faddeev components, 
followed by a second decomposition into Yakubovsky components. For identical bosons, 
one has two Yakubovsky components if) a and if>B- We start from the Yakubovsky equations 
including a general three-body force in the form written down by Glockle and Kamada 



34(. The full four-body bound state wave function is decomposed into the Yakubovsky 



components if) a and ips via 

* = (1 + (1 + P)P M )(l + P)4>A + (1 + P)(l + P)lpB , (11) 

where P^ exchanges particles % and j , P is defined in Eq. (E]) , and P is given by 

P = P13P24 ■ (12) 
The equations for the two wave function components read: 

if) A = G t 12 P[{l + P 34 )^a + M + ^(1 + G t 12 )G V 3 * , 

fa = G t 12 P[(l + P 34 )^A + M , (13) 

where t\2 denotes the two-body t-matrix for particles 1 and 2 and V3 is the three-body 
force defined in Eq. (jSJ). Note that the three-body force couples to the full four-body wave 
function The factor of one third in front of the three-body force term arises because we 
insert the full three-body interaction for V3. This is possible since we consider three-body 
contact interactions which are symmetric under the exchange of any pair of particles. 

The renormalization analysis of the four-body system is complicated by the cutoff depen- 
dence of the number of bound states in the three-body subsystems. The further the cutoff 
A is increased, the more three-body bound states appear. While this has no influence on 
low-energy three-body observables, it creates an instability in the four-body system which 
can collapse into a deep three-body bound state plus another particle. This limits cutoff 
variations to an interval Ao < A < 22.7 Ao for some Ao, in which the number of three-body 
bound states remains constant. Despite this restriction, the cutoff can be varied by more 
than a factor of ten which is sufficient to study the renormalization properties and obtain 
converged numerical results 

In Ref. [23], it was pointed out that the calculation of [I2J for the 4 He tetramer was 
limited to repulsive three-body forces A 3 > 0. The ability to renormalize without a four- 
body parameter was attributed to a very small probability of the four particles to be close 
together due to the repulsive three-body force. We note that the sign of the three-body force 
is cutoff dependent because of the limit cycle behavior of A3 (A) in Eq. fflOl) . In Ref. [jjj], 
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we have performed calculations with both attractive and repulsive three-body forces. In 
particular, the three-body force is attractive for cutoffs close to Ao- However, in this case 
A ~ (m^f) 1 / 2 and the ground state energy is far from the converged value. The final 
converged value for B^ was indeed obtained with a repulsive three-body force. The excited 
state energy B^ , however, is already converged at smaller cutoffs where the three-body force 
is attractive. Increasing the cutoff beyond 22.7 Ao, where the three-body force is attractive 
again would require to project out the unphysical deep three-body state that appears in this 
case. Such a subtraction is beyond the scope of this work. In the remainder of the paper, 
we will follow l3] and not include a four-body parameter in our calculations. 



III. FOUR-BODY UNIVERSALITY 

We now apply this effective theory to calculate the universal properties of the 4-boson 
bound state spectrum and universal scaling functions. We do not discuss the 2- and 3-boson 
spectrum in detail, since such a discussion can be found in Refs. 0, Q, fl3 |. 



A. Bound State Spectrum 

The four-body spectrum is most conveniently discussed in a generalized Efimov plot. 
This plot was introduced by Vitaly Efimov to summarize the universal properties of the 
three-body spectrum 0, 0]. The set of all possible low-energy three-body states can be 
represented as points (a -1 , K) on the plane whose horizontal axis is 1/a and whose vertical 
axis is the momentum variable 

K = sign(E) y/m\E\ . (14) 

It is convenient to introduce polar coordinates consisting of a radial variable H and an 
angular variable £ defined by 

1/a = #cos£, K = Hsm£. (15) 

In terms of these polar coordinates, the discrete scaling symmetry in the three-body system 
is simply a rescaling of the radial variable: H — ► exp(n-K / sq) H. 

The a~ x -K plane for three identical bosons is shown in Fig. [TJ We will refer to the bosons 
simply as particles in the following. The possible states are three-particle scattering states 
(PPP), particle-dimer scattering states (PD), and Efimov trimers (T). The threshold for 
scattering states is indicated by the hatched area. The Efimov trimers are represented by 
the heavy lines below the threshold. 1 There are infinitely many branches of Efimov trimers 
close to threshold, but only a few are shown. A given physical system has a specific value 
of the scattering length and can be represented by a vertical line. The infinite scattering 
length limit corresponds to tuning this line to the K axis. For a real physical system, not 
all states will behave as depicted in the figure. With momenta of order 1/7 (or energies of 



The curves for the trimer binding energies in Fig. [T] actually correspond to plotting H 1 ' 4 sin £ versus 
i? 1 / 4 cos^. This effectively reduces the discrete symmetry factor 22.7 down to 22. 7 1 / 4 = 2.2, allowing a 
greater range of a -1 and K to be shown in the Figure. 
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FIG. 1: The a —K plane for the three-body problem. The allowed regions for three-particle 
scattering states and particle-dimer scattering states are labeled PPP and PD, respectively. The 
heavy lines labeled T are two of the infinitely many branches of Efimov states. The cross-hatching 
indicates the threshold for scattering states. 



order l/(m/ 2 )) one is able to probe details of the short-distance mechanism leading to the 
large scattering length and the effective theory does no longer apply. As a consequence, only 
states with Kl <C 1 and systems with l/\a\ <C 1 will show the universal behavior. 

In Fig. [21 we generalize this plot to the four-body system and plot the square root of 
the four-body energy versus the inverse scattering length. The lines representing the trimer 
energies in Fig. [1] now become the scattering thresholds for trimer-particle scattering. The 
vertical dotted line indicates the limit of infinite scattering length. Our results for the 
four-body ground and excited state energies Bf ] and Bp are given by the circles and 
triangles, respectively. Four-body states (tetramers) can only be stable if their energy is 
below all scattering thresholds, otherwise they become resonances and aquire a width from 
the decay into the corresponding scattering states. The lower (upper) solid lines indicate the 
thresholds for scattering of a ground state (excited state) trimer and another particle. For 
positive a, there are also scattering thresholds for scattering of two dimers and a dimer plus 
two particles indicated by the dash-dotted and dashed lines, respectively. These thresholds 
apply to both tetramer states. A given physical system is again characterized by a vertical 
line corresponding to a particular value of the scattering length. Depending on the value of 
the scattering length, the ground state trimer-particle threshold or the dimer-dimer threshold 
can be closest to the tetramer states. In this paper, we focus on the region where the closest 
threshold is given by ground state trimer-particle scattering. If the plot was extended further 
to the right, the dimer-dimer threshold will eventually become the lowest one. 

We have obtained the dependence of the four-body energies on the scattering length a 
for fixed three-body parameter L3 from explicit solutions of the Yakubovsky equations ([131) . 
All quantities are given in units of L3. Note that L3 is only defined up to a rescaling by the 
discrete scaling factor exp(n7r/s ) ~ (22.7) n , so the absolute scale of Fig. [2] is arbitrary up to 
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a>0 • 

i i i i i i i i 

-0.3 0.3 0.6 

1/(L 3 a) 

FIG. 2: The a^ 1 -K plane for the four-body problem. The circles and triangles indicate the four- 
body ground and excited state energies and while the lower (upper) solid lines give the 
thresholds for decay into a ground state (excited state) trimer and a particle. The dash-dotted 
(dashed) lines give the thresholds for decay into two dimers (a dimer and two particles). The 
vertical dotted line indicates infinite scattering length. All quantities are given in units of the 
three-body parameter L3. 



■1.5 



a < 



such factors. The circles and triangles show our results for the four-body ground and excited 
state energies B^ and as a function of the scattering length, respectively. The tetramer 
excited state remains close to the threshold for decay into a trimer and another particle for 
essentially all calculated values of 1/a. It moves away from this threshold near the point 
where the threshold disappears at K = 0. The tetramer ground state is considerably deeper. 
It stays essentially parallel to the excited state but starts to move away for small positive 
scattering lengths. Due to numerical instabilities, we were not able to reach the region where 
the two-dimer threshold becomes relevant. 

For the limits of universal behavior the same restrictions as in the three-body system 
apply. With momenta of order 1/7 (or energies of order l/(ml 2 )) one is able to probe details 
of the short- distance mechanism leading to the large scattering length and the effective 
theory does no longer apply. As a consequence, only states with Kl -C 1 and systems with 
l/\a\ <C 1 will show the universal behavior in Fig. [2j Which states are within the universal 
window then depends on the actual values of a and / in the system under consideration. For 
4 He atoms there is some evidence that both tetramer states are within this range (l2| . 
It is interesting to note that the universal theory always produces two four-body states 



3JJ . The number of four-body states is independent of the number of three-body states. The 
latter can be expected since the excited three-body states are shallower by at least a factor 
of 515 and therefore should have little influence on the four-body states. The discrete scaling 
symmetry then suggests that there are exactly two four-body resonances between any two 
three-body states with the shallower of the four-body states having a slightly lower energy 
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FIG. 3: The correlation between the three-body ground state energy and the four-body ground 
and excited state energies s| (circles) and (triangles) in units of l/(ma 2 ) for a > 0. The 
inset shows the threshold region in more detail. 

than the shallower of the two three-body states. These states are stable for some range of 
cutoffs, but become resonances with the appearance of a deeper three-body state when the 
cutoff is increased beyond the critical value 22.7Ao discussed in section Hi CI Checking this 
conjecture would require finding resonance poles in the complex plane which is beyond the 
scope of this work. 



B. Universal Correlations 

We now consider the universal correlations in the three-body system for negative scat- 
tering length. In particular, we focus on the correlation between the binding energies in 
the three- and four-body systems. Such correlations were first observed in nuclear physics 
and are known as Tjon lines j35|. For recent discussions of the nuclear Tjon line in the con- 
text of the EFT for short-range interactions and low-momentum nucleon-nucleon potentials, 
sec Refs. [TsUll. The Tjon lines for spinless bosons with positive scattering length were 



discussed in jl2[ for the range of binding energies relevant to 4 He atoms. 

Here, we discuss the correlation between the three-body ground state energy Bf and the 
four-body ground and excited state energies and £4 for values of B^ from threshold 
up to B% pa 200/ (ma 2 ). In Fig. [31 we show the correlation between the ground state energy 
in the three-body system and the ground and excited state energies in the four-body system 
for a > 0. The corresponding plot for negative scattering length is shown in Fig. HI 

The inset in both figures shows the threshold region in more detail. The nonlinear 
behavior of the Tjon lines in this region is evident. For positive scattering length, we can 
not calculate all the way down to the three-body threshold at B% = I /(ma 2 ) due to numerical 
instabilities. The correlations have positive curvature for a > and negative curvature for 
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FIG. 4: The correlation between the three-body ground state energy and the four-body ground 
and excited state energies s| (circles) and (triangles) in units of l/(ma 2 ) for a < 0. The 
inset shows the threshold region in more detail. 

a < 0. Outside the threshold region, the correlation is approximately linear independent of 
the sign of a. A similar linear relation holds for the correlation between different three-body 
energies (e.g. ground and excited states) [jjj]. This is a direct consequence of the discrete 
scaling symmetry in the three-body system. For 1/a = 0, this symmetry ensures that the 
binding energies of subsequent three-body states differ by factors of approximately 515. If 
the scattering length is finite, the correlation is still linear to a good approximation. Only 
for states close to threshold the linearity breaks down. In this case, there is still an exact 
scaling symmetry, but it relates states corresponding to different values of the scattering 
length. Since no new parameter enters at leading order in the four-body system, the above 
arguments immediately generalize to the four-body case. Away from the thresholds, the 
binding energies of the ground and excited four-body states are related to the deepest three- 
body state by factors of 5 and 1.01, respectively. 

Similar linear correlations between the binding energies of N- and (N — l)-body systems 
with positive scattering length and N > A were recently found by Hannah and Blume [36]. 
They studied bosonic clusters with up to iV = 40 atoms interacting additively through two- 
body van der Waals potentials. Using exact Monte Carlo methods, they determined the 
energies as well as the interparticle distances. They found approximately linear relations for 
the correlation between the ground state energy ratios B$ / B$_ 2 and Bf_jBf_ 2 for N 
up to 10. This observation suggests that the discrete scaling symmetry holds up to iV = 10. 
Consequently, no new parameters would enter into the effective theory at leading order for 
up to 10 particles. If true, it would open up the exciting possibility to predict the universal 
properties of TV-body system up to N = 10 from the scattering length a and the three-body 
parameter L 3 alone. 



11 



IV. SUMMARY AND OUTLOOK 



In this paper, we have studied the universal properties of the four-body system with 
large scattering length. These properties are useful for understanding the cluster structure 
of nuclei 24j as part of the planned experimental program at FAIR and experiments on 
the creation of weakly-bound tetramers with cold atoms close to a Feshbach resonance. We 
have concentrated on the bound state problem of four bosons starting from the Yakubovsky 
equations. We have constructed an effective interaction potential including both a two- and 
three-body contact interaction. This is the minimal set of contact interactions required for 
renormalization of the three-body problem. The two parameters of the effective potential 
were determined from matching to the binding energy of the dimer (a > 0) or the two-body 
scattering length (a < 0) and the excited state of the trimer. We have then solved the 
four-body bound state problem and obtained the scattering length dependence of the four- 
body bound state spectrum. This dependence was presented in a generalized Efimov plot. 
For all considered values of a, we have found two four-body states. The tetramer excited 
state stays close to the trimer-particle threshold, while the ground state is considerably 
deeper. We have conjectured, that there are always two four-body resonances between any 
two three-body states. As the cutoff is increased, these states come to life as stable bound 
states, but turn into resonances as additional deep three-body states appear when the cutoff 
is increased beyond a critical value. 

An important consequence of the large scattering length is the emergence of universal 
scaling functions. Since there are only two parameters at leading order, a and L 3 , few- 
body observables normalized to the scattering length a must be correlated and follow a 
line parameterized by L 3 . We have calculated the universal scaling functions relating the 
tetramer energies to the trimer ground state energy. A similar correlation between the 
triton and alpha particle energies is known from nuclear physics as the Tjon line [15|, 



Close to the three-body thresholds, these correlations have positive curvature for a > 
and negative curvature for a < 0. Away from the thresholds, they are linear to very high 
accuracy. This linearity can be understood from the discrete scaling symmetry. Similar linear 
correlations between the binding energies of Y- and (Y — l)-body systems with N > 4 were 
recently found from exact Monte Carlo calculations using two-body van der Waals potentials 
(3^ . This observation suggests the absence of Y-body parameters for Y > 4, which would 
allow to predict Y-body observables from two- and three-body input alone. 

The limits of universal behavior in real systems are determined by the values of a and 
I in those systems. With momenta of order 1/7 (or energies of order l/(ml 2 )) one is able 
to probe details of the short-distance mechanism leading to the large scattering length and 
our effective theory does no longer apply. As a consequence, only states with Kl <C 1 and 
systems with l/\a\ <C 1 will show the universal behavior discussed in this paper. Corrections 
are suppressed by powers of l/\a\ and kI where k is the typical momentum scale in the 



process considered (e.g. the binding momentum for bound states). In Ref. [12j, we have 
shown that our results can be applied to the 4 He tetramer. We are not aware of any data 
on four-body systems with large negative scattering length. 

We also stress that the universality of our results relies on the validity of the renormal- 
ization group analysis in Ref. [13]. In this work, renormalization of the four-body equations 
was achieved without the introduction of a four-body interaction. If this result would not 
hold, the properties of the four-bod y sy stem with large scattering length would also depend 
on a four-body parameter (cf. Ref. [23j|). 
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There are a number of directions that should be pursued in future work. First, it would 
be interesting to extend the calculation to smaller positive values of a, where the relevant 
decay threshold for the tetramer states is the dimer-dimer threshold. This would complete 
the generalized Efimov plot. Our conjecture of two four-body resonances between any two 
trimer states could be tested by extending the calculation to complex values of the energy. 
The widths of such states are currently unknown. Moreover, the application of our results 
to cluster structures in nuclei requires the inclusion of the long-range Coulomb force. Work 
in this direction is in progress. The general power counting for four-body forces is still not 
understood. At which order does the leading four-body interaction enter? In the three-body 
system, e.g., the first and second order correction are due to the two-body effective range 



To| . If a similar situation holds in the four-body system, it would be possible to predict low- 
energy four-body observables up to corrections of order (7/|a|) 3 from two- and three-body 
information alone. The extension of the effective theory to calculate four-body scattering 
observables would be very valuable. The knowledge of the dimer-dimer scattering length, 
for example, is important for experiments with ultracold atoms. For the simpler problem of 
fermions with two spin states (where the three-body paramete r does not contribute), the 
dimer-dimer scattering length was recently calculated exactly [371] and using a perturbative 



e-expansion around d = 4 [38 
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